Relative Ages of Globular Clusters 



Thomas H. Puzia 1 

Sternwarte der Universitat Miinchen, Scheinerstrasse 1, D-81679 Miinchen, Germany 



Abstract. Ages of extragalactic globular clusters can provide valuable insights into 
the formation and evolution of galaxies. In this contribution the photometric methods 
of age dating old globular cluster systems are summarised. The spectroscopic approach 
is reviewed with an emphasis of the right choice of age diagnostics. We present a new 
method of quantifying the relatively best age-sensitive spectroscopic index given the 
quality of a data set and a certain theoretical stellar synthesis model. The relatively 
best diagnostic plot is constructed from the set of Lick indices and used to age date 
globular clusters in several early-type galaxies which are part of a large spectroscopic 
survey of extragalactic globular cluster systems. We find that, independently of host 
galaxy, metal-poor ([Fe/H]< —0.8) globular clusters appear to be old (t > 8 Gyr) and 
coeval. Metal-rich clusters show a wide range of ages from ~ 15 down to a few Gyr. 

1 Introduction 

The foundation of understanding galaxy formation and evolution is the access to 
accurate ages of stellar populations which build up galaxies. All galaxies are ex- 
pected to host multiple stellar populations with different ages, but it is uncertain 
when and how they formed. In general, the ages of galaxies are determined from 
photometric and/or spectroscopic parameters which are more or less correlated 
with the mean lifetime of the observed stellar system. Looking at the diffuse 
light of galaxies one observes the luminosity-weighted sum of all contributions 
emitted by a mix of stellar populations with various ages. Hence, it is meaning- 
less to assign a single age to galaxies and impossible to resolve the age spread of 
the underlying stellar populations by looking at the diffuse light alone. 

In the recent past, globular clusters increasingly played the role of accurate 
age tracers. Globular clusters are simple stellar populations (SSPs, containing 
stars of same age and same metallicity), are formed during major star forma- 
tion events, and are easy to observe out to large distances. The easy age dating 
of these witnesses of galaxy evolution makes it possible to assign single star- 
burst events to globular cluster sub-populations and thus reconstruct the star 
formation histories of galaxies. However, there are still fundamental limits of 
accurate age dating: (1) photometric and spectroscopic age diagnostics suffer 
from the well-known age-metallicity degeneracy and (2) stochastic fluctuations 
in bright stellar evolutionary phases in less-massive clusters (< 10 5 Af©) limit 
the accuracy of age determinations. The remedy for the former point is a smart 
choice of photometric colours and/or spectroscopic indices to reduce the age- 
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Fig. 1. Age dating of globular clusters in NGC 4472 taken from |7J. The dots indicate 
the colours and turn-over magnitudes of the two major globular cluster sub-populations 
which are compared to the SSP model predictions. Different CMD and SSP models 
variations were analysed and gave very similar results. The two globular cluster sub- 
populations were found to differ mainly in metallicity. For details see [ffl] 



mctallicity degeneracy. The latter problem can be by-passed by averaging over 
a large sample of globular clusters. 

In the following, the two fundamental techniques of age dating of old glob- 
ular clusters are discussed. We will focus on globular cluster systems which are 
older than 1 Gyr and hosted in early-type galaxies. Finally, first results of our 
spectroscopic survey of extragalactic globular cluster systems are presented. 



2 Photometric Age Dating 

Photometric ages of extragalactic globular clusters are obtained from the com- 
parison of observed parameters with theoretical model predictions. Such models 
are calibrated on Local Group globular clusters for which the ages are known 
from other techniques. However, the globular clusters in the Local Group are not 
always representative for the variety of clusters found outside the Local Group, 
e.g. in early-type galaxies. To obtain meaningful results for clusters with high 
metallicities which cannot be observed locally one uses extrapolations to derive 
their ages. Absolute age predictions vary among SSP models as they depend 
on the ingredients of the model and the choice of local calibrators. It is essen- 
tial to use more than one model when accurate absolute ages are required to 
secure oneself from model-to-model variations. However, relative ages can be 
reliably obtained using any SSP model as the differential predictions are less 
affected by the choice of a specific model. Photometric age dating methods use 
colour-magnitude diagrams (CMDs) and colour-colour diagrams as their tools. 




Fig. 2. (V - I) vs. (V - K) diagnostic plots for globular clusters in NGC 3115 {left 
panel) and NGC 4365 {right panel) taken from [9. A significant intermediate-age 
globular cluster sub-population is detected in NGC 4365. All used SSP models indicate 
that these objects have age between 2 and 8 Gyr and metallicities ~ O.5Z0 — 3Zq. See 
contribution of Hempel et al. (this volume) and j^] for further systems 

2.1 Colour-Magnitude Method 

The first method is based on mean colours and turn-over magnitudes of major 
globular cluster sub-populations. The requirement is that the sub-populations 
can be separated by their colours and their luminosity function turn-over mag- 
nitudes can be observed. In other words, this method is only applicable to old 
globular cluster systems with "evolved" globular cluster luminosity functions. 
An application is shown in Figure [l] where relative ages of the two major globu- 
lar cluster sub-populations in NGC 4472 can be read off directly from the model 
grid. The two sub-populations were found to be coeval within the errors and 
to differ mostly in metallicity. The use of other photometric filters can enhance 
the predictive power of this method (see contribution of Andres Jordan in this 
volume). The maximum age resolution of this method at old ages (> 10 Gyr) 
can reach 2 — 3 Gyr and is strongly dependent on the photometric quality. This 
technique, however, is premised on the assumption that all sub-populations fea- 
ture the same globular cluster mass function. On the other hand, if extended to 
a wider wavelength range this technique can be used to study variations in the 
globular cluster mass function. For example, accurate V and I photometry in 
NGC 4472 already allowed to constrain the variations in the mass characterizing 
the turn-over to be less than 20%. 

2.2 Colour-Colour Method 

Colour-colour diagrams offer an alternative approach to photometric ages with- 
out the assumption of same underlying globular cluster mass functions. Com- 
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bincd with a wide wavelength coverage from the optical to the near-infrared, 
our group uses the two-colour diagram technique to derive relative ages of glob- 
ular clusters in early-type galaxies (see contribution of Maren Hempel). The 
advantage of linking optical and near-infrared colours is that it reduces the age- 
metallicity degeneracy very efficiently. The reason for this is the enhanced metal- 
licity sensitivity of optical/near- infrared colours (e.g. V — K) due to the sam- 
pling of metallicity-sensitive giant branch luminosity. At similar age-sensitivity as 
purely optical colours (e.g. V — I) this combination minimizes the age-metallicity 
degeneracy and increases the total age-sensitivity in a colour-colour diagram. 

Figure || shows the results of an optical/near-infrared study of two nearby 
early-type galaxies, NGC 3115 and NGC 4365. Our group detected a significant 
population of intermediate-age globular clusters in NGC 4365 || which could 
not be detected from optical colours alone. The results were recently confirmed 
H with spectroscopic methods discussed below. Counterparts were also found in 
NGC 1316 (| and NGC 5846 (see M. Hempel et al. in this volume). The subtle 
but important point of this study is the fact that an intermediate-age stellar 
population can be well-hidden in the diffuse light of galaxies. No hint for an 
intermediate-age stellar population in the diffuse light of NGC 4365 was found 
so far H . We clearly require more near- infrared data on globular cluster systems 
to quantify the mass fractions of such intermediate-age globular cluster sub- 
populations. Future work Q could provide such numbers and strictly constrain 
the importance of certain galaxy formation scenarios. 

3 Spectroscopic Age Dating 

As in the photometric approach, spectroscopic ages are derived from the compar- 
ison of observed parameters with theoretical model predictions. A commonly- 
used tool is the so-called diagnostic plot in which a spectroscopic metallicity 
tracer is plotted against a spectroscopic age indicator. The best spectroscopic 
age indicators for old stellar populations are Balmer lines. However, there is still 
no ideal set or combination of spectroscopic parameters which are entirely free 
from the age-metallicity degeneracy. Yet, combined with a metallicity indicator 
which is almost independent of age robust estimates can be achieved. One way 
to measure the strength of Balmer lines in distant stellar populations is the use 
of spectroscopic indices which do not require high-resolution and high-S/N spec- 
tra. The Lick system jll] defines a number of such indices which are designed 
to measure the strength of three Balmer lines and many metal lines. 

3.1 The relatively best Metallicity Indicator 

Depending on the star formation history the chemical composition of a stellar 
population is subject to abundance variations, in particular the [a/Fe] ratio 
[jl2| . These variations blurr the accuracy of any metallicity diagnostic. New 
SSP models for stellar populations with constant [a/Fe] ratio |l(| allow to 
account for varying abundance ratios and to quantify the effect of changing 
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Fig. 3. Two iso-age lines (3 and 13 Gyr) are plotted as a function metallicity (—2.25 < 
[Z/H] < +0.67) parameterized for three [a/Fe] ratios (0.0, +0.3, and +0.5 dex) in three 
different diagnostic plots, (a): H/3 vs. Mgb. The sensitivity to [a/Fe] can change the 
age prediction by up to ~ 5 Gyr. (b): H/3 vs. (Fe). Sensitivity to [a/Fe] as in panel a. 
(c): H/3 vs. [MgFe]. The isochrones in this diagnostic plot are least affected by [a/Fe] 
variations. Other Balmer line indices show similar behavior. The models are taken from 

m 

[a/Fe] on spectroscopic indices. A good metallicity indicator can now be con- 
structed with the requirement of being insensitive to abundance variations, 
but simultaneously trace the total metallicity. Gonzalez || found discrepant 
age predictions of theoretical models depending on the choice of Lick indices 
used as metallicity indicators. He speculated that this effect might be due to 
variation in the [a/Fe] ratio and could be cured by using the composite in- 
dex [MgFe]= y/Mgb ■ (Fe5270 + Fe5335)/2. In Figure | the effects of changing 
[a/Fe] ratios on old and young isochrones are shown. Indeed, the composite in- 
dex [MgFe] is only little affected by abundance variations and should be used a 
metallicity indicator in all spectroscopic diagnostic plots. 

3.2 The relatively best Age Indicator 

The Lick system defines five Balmer line indices, H/3, H'Ja, ^If, H$a, and HSf- 
The age sensitivity is different for each index and has never been quantified. 
Various combinations of Balmer indices and metallicity indicators are used in 
the literature and assigned equal importance. Here we define a quantity which 
can be used to determine the relatively best Lick Balmer-line index to age date 
stellar populations given the quality of data and any SSP model. The relative 
age-sensitivity of Balmer indices is a function of: 

• rj: mean error of the data 

• £: transformation accuracy to the Lick system 

• 7: mean error of the original Lick spectra 

• 5: accuracy of the Lick fitting functions |U5), O] 

• 2?z : dynamic index range at a given metallicity 

• iSz jt : sensitivity to age and metallicity at a given metallicity and age (i.e. the 
impact of the age-metallicity generacy) 
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Table 1. Summary of the coefficients in equation |2j Column f gives the mean data 
quality of the sample in Figure ||. The coefficients in columns 2-8 are given in units 
of A. Columns 9-12 are given in dex/Gyr while the unit of 1Z in the last column is 
A-dex/Gyr. The average dynamic age range X>z is evaluated at [Fe/H]= —1.35 and 0.0 
between the 1 and 15 Gyr isochrone. S is the mean of the degeneracies at two different 
metallicities [Fe/H]= —1.35 and 0.0 at two different ages 3 and 13 Gyrs. 
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where 5z,t, the degeneracy parameter, is 



Sz.tCO = 



Of 
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We define the quantity 
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• Si ; 



yV + C 2 + 7 2 + <5 2 



(2) 



which is maximal for the best age-sensitive index and is essentially the dynamic 
scale of an index / at a given age and metallicity expressed in units of the 
total uncertainty. Mean values of the dynamic range (2?) and the degeneracy 
parameter (S) are used to evaluate 1Z. Both former parameters are derived from 
the SSP models of and summarised in Table |l|. 

We find that from our data the relatively best age diagnostic is the H7A index 
followed by the indices H^, H/3, Hyp, and H5f- If our data set were infinitely 
accurate (i.e. 77 = in Eq. |^) the order of the Balmer-indices would remain the 
same. However, 1Z is subject to change when using different SSP models. 



3.3 Spectroscopic Ages of Globular Clusters in Early-type Galaxies 

First results of the age dating of globular clusters from our spectroscopic survey 
of globular cluster systems in early-type galaxies are shown in Figure I Inde- 
pendent of galaxy, metal-poor ([Fe/H]< —0.8) globular clusters appear to be 
old and coeval. Metal-rich ([Fe/H]> —0.8) globular clusters, on the other hand, 
cover a larger age range and appear on average younger. This interesting re- 
sult indicates that metal-poor globular clusters might have preferentially formed 
very early, perhaps during or even before the assembly of their host galaxies. 
The metal-rich sub-population appears to form continuously until the recent 
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Fig. 4. Globular clusters in the early-type galaxies NGC 1380, NGC 2434, NGC 3115, 
and NGC 3379. SSP models were taken from Isochrones (solid lines) are plotted 
for the ages 1 to 15 Gyr. A very young iso-age line for a stellar population of 0.8 Gyr 
is also plotted (dashed lines). Iso-metallicity lines (dotted lines) are plotted for values 
-2.25 < [Z/H] < +0.67 

past. This in turn implies that only the metal- rich population of clusters is aug- 
mented in later star formation events (see also M. Kissler-Patig et al. in this 
volume and ||). 

Some caveats have to be addressed when ages are derived from spectroscopic 
diagnostic plots. In the metal-poor regime ([Fe/H] < —0.8) at old ages (t > 10 
Gyr), SSP model tracks start to overlap due to the increasing importance of hot 
blue horizontal branch stars. The consequence are ambiguous age predictions 
for old metal-poor populations and a potential age spread at old ages which 
cannot be ruled out. In the metal-rich regime ([Fe/H] > —0.8) the accurate age 
determination might also be hampered by an unknown contribution of the blue 
horizontal branch as a consequence of inaccurate modelling of the mass loss on 
the giant branch. 
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